The Difficulty of Gate Control in Molecular Transistors 
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The electrostatic gating effects on molecular transistors are investigated using the density func- 
tional theory (DFT) combined with the nonequilibrium Green's function (NEGF) method. When 
molecular energy levels are away from the Fermi energy they can be linearly shifted by the gate 
voltage, which is consistent with recent experimental observations [Nature 462, 1039 (2009)]. How- 
ever, when they move near to the Fermi energy (turn-on process), the shifts become extremely small 
and almost independent of the gate voltage. The fact that the conductance may be beyond the gate 
control in the "ON" state will challenge the implementation of molecular transistors. 



Introduction. — The molecular transistor is the most 
essential element of molecular electronics. How to control 
molecular orbitals by external gate voltage {Vq) is the key 
problem of molecular transistors. Although two-terminal 
molecular devices (with the source and drain electrodes) 
have been extensively studied for many years, the real 
three-terminal ones (with the gate electrode added in) 
have only been developed recently in the laboratory by 
Song et al [ij. They found that the transport current 
of Au/l,4-benzenedithiol (BDT)/Au (as well as Au/1,8- 
octanedithiol/Au) in a transistor structure can be di- 
rectly controlled by gate modulating molecular orbitals. 
A linear relationship between the gate voltage and the 
shift of the molecular orbital energy has been also ob- 
served [ij. Based on their experiment, it seems that the 
practical molecular field-effect-transistor (FET) will be 
realized soon. 

From a practical point of view, in order to control the 
carriers or the "conductive channel" by the gate volt- 
age, the molecular FET should work well on three modes 
of operation (as the traditional counterpart does): cutoff, 
linear and saturation regions. We comment that the gate 
control of molecular orbital only observed in the cutoff 
region, i.e. the channel between the source and drain 
keeping in the "OFF" state in the experiment [1]. Of 
course, one expects molecular-orbital-gating occurring in 
the cutoff region is also valid in the other two regions (the 
"ON" states). However, such natural extrapolation in mi- 
cro electrons is questionable in molecular transistors, be- 
cause at the molecular level the motion of carriers obeys 
quantum mechanical which may make the real situation 
completely different. In this letter, we will address this 
issue via the first principles calculations-the density func- 
tional theory (DFT) combined with the nonequilibrium 
Green's function (NEGF) method. We will demonstrate 
that above extrapolation does be problematic: once the 
channel between the source and drain being turned on, 
the conductance of the molecular transistor may be be- 
yond the control of gate voltage. This may be a serious 
challenge to the implementation of molecular transistors. 



Method. — The gate-controlled Au/BDT/Au struc- 
ture is the object of our theoretical study. However, the 
aim of the present work is not to quantitatively simulate 
the experment of Song et al , but to qualitatively reveal 
the physics behind it and the physics being missed in it. 
Based on this point, we apply the gate voltage via an ad- 
ditional external potential instead of adding a real gate 
electrode in the first principle calculations (the latter is 
very difficult for DFT+NEGF calculations). 

The Au/BDT/Au structure consists three parts: the 
left and right electrodes plus the scattering region (with 
gate voltage being applied). Both of the electrodes are 
periodically repeated structures of the super cell formed 
by 3 layers of 3x3 Au(lll) slabs, while the scattering re- 
gion is composed of a single molecule BDT and four layers 
of 3x3 Au(lll) slabs on each side. The BDT molecule 
bounds to the electrode surfaces through the thiol groups 
contacting with Au adatoms on each side, as shown in 
FigurelH It is reported that when BDT molecule contacts 
with Au leads, the structure with the thiol group S-H 
non-dissociated is energetically more favorite P, 111, thus 
we adopted this kind of structure in our calculation. As 
in most literatures, BDT molecule is set perpendicular to 
the Au surfaces. Although the molecule structure and its 
coupling to the electrodes may change the conductance 
of BDT single molecular junction, they should not cause 
qualitatively difference to the results reported here. We 
will neglect those effects and focus on the gate-controlled 
conductance of single-structure molecular junctions. 

All of the calculations are performed with spin- 
polarized DFT implemented in the Siesta code using 
Perdew-Burke-Ernzerhof exchange-correlation functional 
fill. Double-C plus polarization function(DZP) basis set 
is used, together with a mesh cutoff of 200 Ry and norm- 
conserving pseudopotentials. The distance between the 
left and right lead surfaces is relaxed to the equilibrium 
position. Then the BDT molecule together with the near- 
est layer of Au atoms are fully relaxed to fulfill the energy 
and force convergences of 10~^ eV and 0.01 eV/ A respec- 
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Figure 1: (color online) Upper panel: structure geometry 
of Au/l,4-benzeneditliiol (BDT)/Au junctions (the scatter- 
ing region). Lower panel: z dependent gate voltage profile in 
the form shown in Eq.(p^. 

The transport property is calculated using 
DFT+NEGF method implemented in the recent 
Siesta code (the TranSiesta part 0)- The periodic 
boundary conditions are used in directions orthogonal 
to the transport direction. A 3x3 uniform k meshes are 
sampled over this two-dimensional Brillouin zone. The 
density is integrated with 32 points on the semicircle, 
32 points parallel to the real axis, and 16 poles. The 
transmission spectrum at zero bias is obtained by 
integrating over 1000 points. 

To study the electrostatic gating effects, we follow the 
idea of Morari et al.^ i.e. the gate voltage is simulated 
by adding external potential determined solely by the 
coordinate z along the transport direction [8]. In our 
work, the external potential is in the form of 

U{z) = eVG{ ^^J_^^^^^^ - l + e(/-..)/AJ (1) 

where Az = 0.5 A. The parameters zi and Z2 are used to 
ensure that the applied gate-potential lies in the molecule 
region and decrease to zero smoothly when reaching the 
electrodes, as the potential profile shown in Figure [TJ 
Comparing with experiments, this simulation of gate 
voltage obviously lacks the effect of the gate dielectric 
(AI2O3 layer in ReflU). This effect should not distinctly 
change the physics presented in this work, at least in the 
qualitatively level. 

Results and discussions . — We first check our calcu- 
lation before the gate voltage being applied and the re- 
sults are summarized in Figure [21 where the upper panel 
depicts the zero-bias transmission spectrum T{E) (the 
Fermi energy is set as the zero point of energy), with 
its insert showing the calculated I — V curve. One can 
see that the current is of the same order (/iA) with ex- 
perimental measurement [1], which is apparently an ac- 
ceptable result for static DFT. It further confirms the 



chosen junction geometry (see Figure [T]) is reasonable. 
The asymmetric character of / — V curve is also con- 
sistent with experiment [1], which may result from the 
asymmetric coupling of BDT molecule with the left and 
right electrodes. 

From Figure O we find five transmission peaks within 
the energy range of —4 to 4 eV and they are labeled 
as peaks 1-5, where peaks 1-3 (centered on —3.7, —2.9 
and —1.9 eV respectively) correspond to hole-transport 
while peaks 4-5 (centered on 1.2 and 3.0 eV respectively) 
to electron-transport. Actually, we can further identify 
which element (C, S or Au) contributes most to those 
peaks via the spatial charge density [also called the local 
density of states (LDOS)] within energy ranges around 
them, as shown in the lower panel of Figure [2l It can 
be seen that C atoms contribute most to both peaks 1 
and 4, while Au atoms (including those linking BDT and 
bulk electrodes) do most to peaks 2 and 3. The main 
contribution to peak 5 is clearly from S atoms and the 
S-C hybridization. From those information, we can con- 
clude that peaks 1 and 4 are respectively result from 
the conduction of the highest occupied molecular or- 
bital (HOMO) and lowest unoccupied molecular orbital 
(LUMO) of BDT. Since the LUMO is nearest to E^, it 
will dominate the transport properties of the device un- 
der low bias voltage. This result seems not consistent 
with the experiment of Song et al^ where the HOMO is 
observed to be turned on prior to the LUMO. It indi- 
cates that some factors may be missed in our ideal DFT 
calculations, e.g. many body corrections jH. We argue 
that the physics of activating and controlling LUMO or 
HOMO should be same for molecular transistors [11]. 

Now, let us show the complete story when gate voltage 
being applied in the form of Eq.([T]) (cf. Figure [T]). In Fig- 
urelD^a), we plot the transmission spectrum under various 
gate voltages (from— 8 V to +20 V with the interval of 
2 V). As we have demonstrated in Figure O the LUMO 
of BDT plays the leading role in the molecular-orbital- 
gating process under low bias. It should be enough for 
us to focus on the gate control of this orbital, e.g. the 
change of transmission peak 4 with Vg • To be more spe- 
cific. Figure [31(b) depicts the projected density of states 
(PDOS) only on carbon atoms (to highlight the change 
of LUMO) under synchronously changing Vq as shown 
in Figure [31(a). Combining Figure [31(a) and [31(b), we can 
clearly see the modulation of Vq on the energy of the 
LUMO (£'l) and accurately understand its mechanism. 

At first glance. Figure [3l does confirm the fact that 
the conductance and molecular orbital of BDT junction 
can be directly controlled by Vq- However, in different 
range of Vq^ molecular-orbital-gating exhibits different 
behavior patterns. Generally speaking. Figure [3l [both 
(a) and (b)] can be divided into three regions in term of 
Vg, which are explained in details as follows. 

Region I (— 8 V <Vg < 3 V): In this region, the molec- 
ular transistor is in the "OFF" state due to the near zero 
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Figure 2: (color online) (a) The transmission spectrum T{E) 
for BDT junction under zero gate voltage in the unit of quan- 
tum conductance Go(= 2e^/h) . The insert is the calculated 
I — V curve, (b) The spatial charge density within energy 
ranges around four transmission peaks labeled by 1-5 in (a). 

DOS and near zero T{E) at E^. Actually, it is the re- 
gion where Song et al observed molecular-orbital-gating 
in their experiment [1]. Our theoretical results verify 
their observation that the gate control of molecular or- 
bital is linear and effective. To demonstrate this point 
in more details, we extract the change of the energy of 
the LUMO (defined as Ei, — E-p) as functions of Vq from 
Figure [3] (in all range of Vq^ not only in Region I), and 
summarize the results in Figure IH In that figure, the re- 
markable linear-modulation of Vq on Ei, is clearly shown 
in region I, with the slope being defined as the gate effi- 
ciency factor, a = AEi^/eAVc (cf. Ref,l). From Figure 
m we work out a ~ —0.32 in Region L 

By referring [31(b), one can find a similar linear-control 
of Vg on the energy of the HOMO (^h) in the energy 
range of £^ < E^. The gate efficiency factor, an (= 
AEn/eAVc) ~ 0.32, is the same as that of gate control 
of El, in spite of a sign difference. That result confirms 
our argument that activating and controlling LUMO or 
HOMO should share the same physics. Not surprisingly, 
the activation of the HOMO is much harder than that of 
the LUMO in our theory. 



Our theoretical \a\ (~ 0.32) is comparable with the 
experimental value a ~ 0.22 pL]- That provides a reli- 
able basis for further understanding theoretical results in 
Region H, which has not been reported in experiments 
yet. 
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Figure 3: (color online) Calculated transmission spectrum 
T{E) (a) and projected density of states (PODS) on carbon 
atoms (b) for BDT junction under various gate voltages. 

Region H (3 V <Vg < 16 V): In this region, the molec- 
ular transistor is in the "ON" state due to the finite DOS 
and finite T{E) at E^. As shown in Figure [3] and HI 
the behavior of molecular-orbital-gating in this region is 
completely different from that in Region L both DOS and 
T{E) near E^ change very slowly with further increasing 
Vq once the LUMO has been activated to be conductive. 
If we still linearly fit the gate modulating on in Figure 
m we will get a near zero a (~ —0.005). It indicates that 
the energy of LUMO becomes beyond the control of gate 
voltage in this region. 

The weak dependence of T{E) on Vq means the source- 
drain current /sD at finite bias also independent on Vq, 
i.e. Vg will lose the control of /sd in molecular transistors 
in the "ON" state. That will challenges the implementa- 
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tion of molecular transistors. As we know, /sd at certain 
source-drain bias (Vsd) is supposed to increase with Vq 
linearly in the linear region and quadratically in the sat- 
uration region. Neither of them can be achieved when 
the molecular orbital being beyond the control of Vq- 

In order to gain an insight into the difficulty of gate 
control in Region II, we show the x — y planar-averaged 
converged Hartree potentials (/7h) under various Vq in 
the insert of Figure IH where the peaks correspond to 
chemical bonds while the valleys to chemical groups (CH 
or SH) or Au atoms. By comparing Figure [Hand Figure 
[H we find two potential nodes at the interface between 
BDT and Au electrodes (S-Au bonds). The potential of 
S atoms inside the nodes drops rapidly with increasing 
Vg, while that of Au atoms outside of the nodes grows 
fast. 

Since negative potential tends to trap electrons but 
positive one trap holes, those two nodes indicate that Vq- 
induced negative and positive charges localized at the in- 
terfaces separates from each other (see Table I for further 
details). As a consequence, an inner electric field builds 
up to offset the Fc-induced one. This process is very 
similar to the shielding of external field by a plate capac- 
itance, we thus called the effective capacitance near the 
(left or right) interface as structural capacitance (Cs)- It 
can explain why the effective potential acting on molec- 
ular energy is less than the applied one (|a| < 1), even 
no gate dielectric structure being included. 

Table I: Gate voltage Vg induced changes of charges at left 
(Aul — S^,) and right interfaces {Sr — Aur) and that within 
the total scattering region. Vg varies from— 8 V to 20 V with 
the interval of 4 V. 



Vg 


-8 -4 4 8 12 16 20 


Sl 


0.09 0.04 -0.03 -0.07 -0.11 -0.16 -0.21 


Aul 


-0.13 -0.06 0.08 0.18 0.27 0.37 0.46 


Sr 


0.11 0.05 -0.06 -0.11 -0.17 -0.23 -0.29 


Aur 


-0.16 -0.07 0.08 0.18 0.29 0.40 0.50 


Total 


0.67 0.32 -0.41 -0.94 -1.48 -2.02 -2.54 



When the molecular transistor is in the "ON" state, the 
finite DOS at [^(£^f)1 will lead to the so called quan- 
tum capacitance effects [9i,[lo|- The quantum capacitance 
Cq can be defined through Cq = c^D^Ey)- In a realistic 
transistor, the effective capacitance C is the series com- 
bination of the electrostatic capacitance {Ce) and Cq. 
In our theory, C^; is replaced by structural capacitance 
Cs^ therefore, C = CsCq / {Cs ^Cq)^ which is obviously 
dominated by the smaller of the two. As shown in the 
insert of Figure HI the narrow distance between oppo- 
site charges indicates that Cs should be very large (plate 
capacitance model), which makes Cq dominate. As a 
consequence, Cq effectively reduces the external poten- 
tial acting on BDT, or even completely shield Vc-induced 
field to make Vq loss control of molecular orbit als. 
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Figure 4: (color online) The change of molecular orbital 
(LUMO of BDT) energy as a function of the gate voltage. 
The straight lines are eye- guide lines with different slopes. 
The insert is the x — y planar- averaged converged Hartree 
potentials under gate voltage increasing from— 8 V to +20 V 
with the interval of 2 V, as indicated by the arrow. The dif- 
ficulty of gate control in molecular transistors in the "ON" 
state is clearly shown in the figure. 



Region III (16 V <Vg < 20 V): The molecular transis- 
tor is turned back to "OFF" state when sufficient large 
Vg shifts the LUMO of BDT away from E^. The efficient 
linear-modulation of Vq on molecular orbitals is also re- 
stored but with a different slope {a ~ —0.07) from that in 
Region I. When the LUMO+1 of BDT is shifted near E^, 
one can expect a similar process as happens in Region II. 
In that case, the electron-electron correlation should be 
carefully treated (beyond the mean-field theory) due to 
more electrons being trapped within the BDT molecule, 
which has been already beyond the scope of the present 
work. 

Summary. — In summary, we have studied the gate 
control of molecular orbitals by DFT+NEGF method. 
When molecular energy levels are away from the Fermi 
energy they can be linearly shifted by the gate voltage, 
which is consistent with recent experimental observa- 
tions. When they move near to the Fermi energy, the 
shifts become extremely small and almost independent 
of the gate voltage, which indicates that it may be diffi- 
cult for the gate voltage to control molecular orbitals in 
the "ON" state. 
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